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NOMENCLATURE 


c 

chord  length 

1 

span  of  the  cascade  blade 

Ly 

coordinate  of  tip  vortex  center  from  wall 

U 

coordinate  of  tip  vortex  center  from  suction  surface 

R 

vortex  core  radius 

Rn 

Reynolds  number  (Vref  c  )  /  v 

s 

pitch  of  the  casade  blade 

V.ef 

magnitude  of  the  velocity  at  inlet  midspan 

Vx 

velocity  component  in  X  direction 

X 

distance  from  the  leading  edge  of  the  blade,  in  X  direction 

X,Y,Z 

Cartesian  coordinate  system 

V 

kinematic  viscosity 

ABBREVIATIONS 

CPU 

Central  Process  Unit 

NACA 

National  Advisory  Committee  for  Aeronautics 

NAVSEA 

Naval  Sea  Systems  Command 

NSWCCD 

Naval  Surface  Warfare  Center,  Carderock  Division 

FANS 

Reynolds  Averaged  Navier-Stokes 

SGI 

Silicon  Graphics  Incoperated 

mm 

millimeter 

m/s 

meter  per  second 
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ABSTRACT 


The  evolution  of  the  tip  separation  vortex,  the  tip  leakage  vortex  and  the 
passage  vortex  generated  by  flow  over  a  low-speed  linear  compressor  cascade 
with  tip  clearance  was  studied  by  solving  the  Navier-Stokes  equations 
numerically.  The  discretization  is  based  on  flux-difference  splitting,  implicit 
high-resolution  scheme.  Up  to  a  third-order  accurate  upwind  differencing  is 
achieved  through  a  reconstruction  of  the  solution  from  its  cell  averages.  Balwin- 
Lomax  turbulence  model  with  blending  function  to  account  for  the  corner  flow 
was  used  for  the  turbulence  closure.  The  formation  and  the  development  of  the 
vortices  were  simulated  and  were  compared  with  experimental  observation. 

ADMINISTRATIVE  INFORMATION 

This  work  was  sponsored  by  the  Advanced  Submarine  R&D  office  (93R),  Naval  Sea 
Systems  Command  (NAVSEA)  and  administrated  by  the  Advanced  Propulsor  Development 
Program  Office  (Code  508),  Naval  Surface  Warfare  Center,  Carderock  Division  (NSWCCD).  The 
work  was  conducted  by  the  Hydromechanics  Directorate,  Propulsion  and  Fluid  Systems 
Department  (Code  5400),  NSWCCD,  imder  work  unit  number  1-5080-172-54. 

INTRODUCTION 

A  structure  of  the  flow  in  the  blade  row  of  a  turbomachine  is  predominately  three- 
dimensional.  Vortices  are  generated  near  the  blade  surface  and  within  the  tip  gap.  The  vortices  of 
particular  interest  are  the  leading  edge  separation  vortex,  the  leakage  vortex  and  the  passage 
vortex.  Excessive  activities  of  vortices  may  cause  some  detrimental  effects,  such  as  efficiency 
loss,  noise  generation,  cavitation  and  blade  erosion.  Numerous  experiments  have  been  conducted 
to  gain  an  imderstanding  of  the  nature  of  the  vortex  flow  [1,2,3,4,5,6,7,8].  Complementarily, 
numerical  methods  have  been  developed  to  facilitate  predictions  [2,9,10,1 1,12],  so  that  an 
optimal  design  can  be  achieved. 

The  flow  visualizations  and  measurements  fi'om  the  tip  gap  experiments  conducted  by 
Kang  and  Hirsch  [3,4,5]  clearly  indicated  that  a  multiple  vortices  structure  existed  around  the 
blade  tip.  From  the  numerical  simulation  point  of  view,  this  experiment  provides  an  ideal  and 
challenging  environment.  Recent  numerical  simulations  of  the  vortex  related  flows  [2,9,13] 
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suggested  that  with  a  dense  grid  distribution  and  a  high-order  differencing  technique  detailed 
flow  structure  within  the  vortex  core  can  be  qualitatively  simulated. 

The  purpose  of  the  present  study  is  to  examine  the  feasibility  of  simulating  the  multiple 
vortices  structure  observed  by  Kang  and  Hirsch  with  a  third-order  accurate  upwind  finite- 
differencing  high-resolution  Reynolds-averaged  Navier-Stokes  (RANS)  solver  [14,15], 

DESCRIPTION  OF  THE  EXPERIMENT 

The  detailed  description  of  the  experimental  facility,  the  setup  and  the  offset  of  the  blade 
section  can  be  found  in  reference  3.  Only  a  brief  and  computationally  relevant  outline  is  given 
here.  The  cascade  consists  of  seven  NACA  65-1810  blades  having  one  free  ends  with  flat  tips. 
The  aspect  ratio  is  1 .0  with  blade  chord  length  of  200  mm.  The  pitch  is  1 80  mm.  The  tip 
clearance  is  2.0  percent  of  the  chord  length  (4mm),  with  an  uncertainty  of  ±0. 1mm  for  the  seven 
tip  gaps.  The  blade  stagger  angle  is  10".  The  blade  inlet  angle  is  32.5"  and  the  blade  outlet  angle  is 
—12.5  .  The  mass  averaged  air  inlet  and  exit  angles,  measured  40  %  of  chord  length  upstream  and 
25  %  of  chord  length  downstream  are  29.3°  and  -2.5"  (design  values  are  30"  and  -4.5°, 
respectively).  The  tests  are  conducted  with  the  inlet  flow  velocity  being  kept  constant  at  23.7  m/s. 
The  test  Reynolds  number  Rn,  based  on  the  inlet  velocity  and  the  blade  chord  length  is  about 
2.9x  10^.  The  free  stream  turbulence  intensity  is  3.4%. 

A  five-hole  probe  (the  diameter  of  the  sensitive  part  of  the  probe  is  2.6  mm)  is  used  to 
acquire  the  three-dimensional  flow  ahead  of,  within  and  behind  the  cascade.  Surface  static 
pressures  (on  both  the  wall  and  blade  surfaces)  are  recorded  by  static  pressure  taps.  The  five-hole 
probe  is  calibrated  in  a  jet.  The  accuracy  of  the  measured  velocity  and  total  pressure  are  1%  of  the 
inlet  mid-span  values  where  all  the  reference  parameters  are  taken. 

NUMERICAL  SCHEME 

The  numerical  solver  is  based  on  the  three-dimensional  incompressible  Navier-Stokes 
equations  expressed  in  the  conservative  law  form.  The  equations  are  coupled  with  the  continuity 
equation  through  the  artificial  compressibility.  The  diffusive  terms  are  discretized  with  a  central 
differencing  scheme,  and  the  convective  terms  are  discretized  with  the  Roe’s  flux-difference 
splitting  scheme  [16].  Third-order  accuracy  is  achieved  by  reconstruction  of  the  solution  from  its 
cell-averages  [14,15].  The  reconstruction  is  accomplished  by  linear  interpolation  of  the 
characteristic  variables  in  the  computational  domain.  The  non-TVD  (Total  Variation 
Diminishing)  forms  are  used  in  order  to  avoid  excessive  diffusion  in  the  vortex  core  region.  Time 
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marching  is  performed  with  first-order  accurate  Euler  backward  time  differencing.  The  Modified 
Baldwin-Lomax  Model  for  three-dimensional  flow  formulated  by  Chima  et  al.  [17]  is  used  for 
turbulence  closure. 


GRID  GENERATION 

The  Cartesian  coordinates  system  is  used  to  describe  the  cascade  geometry.  The  X  axis 
runs  from  upstream  towards  downstream,  the  Z  axis  points  from  the  root  of  the  blades  towards 
the  tip,  and  the  Y  axis  is  chosen  along  the  pitchwise  direction.  The  origin  of  the  coordinates 
system  is  located  on  the  blade  leading  edge  at  the  root.  The  orientation  of  the  blade  is  defined  by 
its  stagger  angle  ( 1 0"). 

The  computational  domain  consists  of  two  portions,  one  part  is  the  passage  bounded 
partially  by  the  pressure  side  surface  of  one  blade  and  the  suction  side  surface  of  the  adjacent 
blade  in  addition  to  the  periodic  boundaries,  and  the  other  part  is  the  complementary  tip  gap 
region.  The  inlet  plane  is  located  three  chord-lengths  upstream  of  the  blade’s  leading  edge,  and 
the  outlet  plane  is  located  three  chord-lengths  downstream  of  the  blade’s  trailing  edge.  Based  on 
the  H-H  type  topology,  grid  systems  for  blade  passage  and  tip  gap  are  generated  with  the  Gridgen 
code  (version  13.4).  In  the  pitchwise  direction,  the  grids  are  clustered  near  the  blade  surfaces.  In 
the  streamwise  direction,  the  grids  are  clustered  around  the  leading  and  the  trailing  edges.  The 
grids  are  distributed  according  to  the  hyper-tangent  rule.  The  first  grid  spacing  off  the  solid 
surface  is  set  to  be  1  x  10'^,  in  order  to  provide  sufficient  resolution  for  the  turbulent  viscous  flow. 
At  the  given  Reynolds  number  (2.9x10^^),  the  first  turbulence  coordinate  y^  off  the  solid  surface 
is  about  1.  Grid  systems  with  three  different  sizes  are  generated  to  represent  the  passage  zone. 
The  sizes  are  195x65x113  (in  axial,  pitchwise,  and  streamwise  directions,  respectively), 
195x107x113,  and  195x121x123.  The  grid  size  for  the  complementary  tip  gap  zone  is 
1 1 1  x33x33.  For  each  of  the  grid  systems,  the  total  grid  nodal  points  are  1 .5  million,  2.5  million 
and  3.0  million,  approximately.  The  increasing  grid  nodal  point  munber  is  to  support  finer 
resolution  of  the  vortex  flow. 


COMPUTATIONAL  RESULTS 

The  results  presented  here  were  obtained  from  the  computations  based  on  the  grid  system 
of  3.0  million  nodal  points.  Computations  are  carried  out  with  6  SGI  Origin  2000  CPUs  in 
parallel.  Figure  1  shows  the  computed  particle  tracing  of  the  vortices  near  the  tip  region.  Figure  2 
shows  the  schematic  of  the  vortex  structure,  reconstructed  through  flow  visualization  from 
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experiment  [4],  Comparison  between  Figures  1  and  2  indicates  that  formation  of  the  observed  tip 
separation  vortex  and  the  tip  leakage  vortex  are  successfully  simulated.  Figure  3  shows  the  ink 
trace  visualization  on  the  blade  tip  surface  [4],  and  Fig.  4  shows  the  corresponding  reconstructed 
schematic  of  the  flow  pattern  [4].  Figure  5  shows  the  computed  limited  skin  friction  lines  on  the 
tip  surface.  Comparison  among  Figures  3-5,  it  can  be  concluded  that  the  positions  of  the 
separation  line  (convergent  line)  and  the  reattachment  line  (divergent  line)  are  predicted  correctly. 
This  is  an  indication  that  the  essential  flow  pattern  at  tip  gap  is  simulated  correctly.  Figure  6 
shows  the  cross  flow  vectors  and  the  particle  traces  in  the  tip  gap  located  at  46%  chord  (x/c  = 
0.46).  Flow  enters  from  the  pressure  side  (left)  and  exits  to  the  suction  side  (right).  Particle 
tracing  shows  that  the  flow  separates  and  then  re-attaches  to  the  tip  surface  near  the  suction  side, 
a  separation  bubble  is  formed.  The  bubble  marks  the  passage  of  the  tip  separation  vortex.  At  the 
leading  edge  of  the  bubble,  two  secondary  separations  can  be  seen.  Along  the  solid  surfaces,  the 
boundary  layers  are  developed.  The  shear  layers  with  embedded  vorticity  are  then  discharged  into 
the  suction  side  of  the  blade.  The  rolling  up  of  the  entrapped  vorticity  constitutes  the  leakage 
vortex.  The  length  of  the  reference  vector  at  the  up  right  of  Figure  6  is  equivalent  to  the 
magnitude  of  the  inlet  velocity  at  mid-span.  Figures  7-10  show  the  cross  flow  velocity  vectors 
and  static  pressure  contours  around  the  tip  leakage  vortex  core  on  the  various  traverse  planes  at 
axial  locations  x/c=0.6,  0.8,  1 .00  and  1 .125,  respectively.  The  length  of  the  reference  vector  at  the 
top  of  the  figures  represents  the  magnitude  of  the  inlet  velocity  at  midspan.  The  location  of  the 
center  of  the  vortex,  where  the  cross  flow  vector  vanishes,  nearly  coincides  with  that  of  the 
minimum  static  pressure.  In  a  previous  grid  refinement  study  conducted  by  Dacles-Mariani  et  al 
[2],  it  was  shown  that  an  average  of  15-20  grid  points  (in  each  of  the  pitchwise  and  the  spanwise 
directions)  in  the  vortex  core  region  are  needed  in  order  to  adequately  resolve  the  flow.  Hsiao  and 
Pauley  drew  similar  conclusion  in  their  study  of  the  tip  vortex  flow  generated  by  a  marine 
propeller  [9].  In  the  present  study,  the  grid  resolution  at  the  core  region,  typically  as  indicated  by 
the  vectors  in  Fig.  8,  is  finer  then  the  20x20  (pitchwise  and  spanwise  directions)  distribution  as 
recommended.  Figure  1 1  shows  the  vorticity  magnitude  contours  at  various  traverse  planes  at 
axial  locations  from  x/c=0.20  to  x/c=1.00  with  an  increment  of  10%  chord  length.  It  shows  that 
the  initial  vortex  roll  up  occurs  around  x/c=0.20,  and  the  vortex  becomes  fully  developed  at 
x/c=0.60.  At  this  point  the  shape  of  the  vortex  becomes  nearly  axisymmetric.  The  vorticity 
magnitude  attends  its  maximum  values  at  x/c=0.8,  and  starts  to  decline  further  downstream. 
Ramaprian  and  Zhang  [6]  suggested  that  as  the  rolling-up  process  continues,  several  competing 
processes  determine  the  peak  vorticity  magnitude  within  the  vortex.  The  first  of  these  is  the 
continuous  trapping,  by  the  vortex,  of  the  vorticity  from  the  shear  layer.  This  tends  to  increase  the 
magnitude  of  the  negative  vorticity  with  distance.  On  the  other  hand,  the  second  process,  namely. 
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diffusion  of  the  vorticity  by  viscosity  and  turbulence,  tends  to  decrease  the  peak  magnitude  of  the 
vorticity  within  the  vortex.  Figure  1 1  also  indicates  that  near  the  trailing  edge,  just  off  the  blade 
surface,  a  vortex  similar  to  the  secondary  vortex  as  shown  in  Fig.  2,  is  developing.  This  secondary 
vortex  can  also  be  traced  from  the  vector  plot  in  Fig.  9.  Figure  12  displays  the  perspective  view  of 
the  tip  vortices  and  the  static  pressure  contours  around  their  cores  at  various  traverse  planes  from 
x/c=0.4  to  x/c=l  .0  with  an  increment  of  10%  chord  length.  The  contours  show  that  at  the  initial 
rolling-up  stage  the  static  pressure  at  the  center  of  the  core  drops;  after  the  vortex  is  fully 
developed,  the  static  pressure  starts  to  recover.  This  recovery  begins  about  at  x/c=0.80,  the  same 
location  as  the  vorticity  magnitude  reaches  its  maximum  values  in  as  shown  in  Fig.  1 1.  Figure  13 
shows  the  cross  flow  streamlines  at  x/c=0.80.  The  left  boundary  is  on  the  suction  side  surface  of  a 
blade  and  the  right  boundary  is  on  the  pressure  side  surface  of  the  next  blade.  Two  vortices  are 
shown  near  the  top.  The  one  near  the  suction  side  surface  is  the  tip  leakage  vortex,  its  shape  is 
nearly  axisymmetric.  The  other  is  the  tip  passage  vortex.  There  are  also  two  vortices  near  the 
root.  One  is  the  comer  vortex  and  the  other  is  the  root  passage  vortex.  Kang  and  Hirsch  illustrated 
such  vortices  at  root  schematically  as  shown  in  Fig.  15  [3].  Detailed  discussion  of  the  formation 
of  these  vortices  can  be  found  in  reference  3.  Figure  14  shows  the  surface  static  pressure 
distribution  on  the  blade  tip  and  on  the  suction  side.  There  is  an  indication  of  pressure  trough  in 
Fig.  14,  near  the  tip  on  the  suction  side  surface.  The  pressure  trough  is  due  to  the  influence  of  the 
low  pressure  in  the  center  of  the  tip  leakage  vortex.  Similar  but  less  intensive  pressure  troughs 
can  be  found  near  the  root,  these  are  caused  by  the  comer  and  root  passage  vortex  as  shown  in 
Fig.  13.  Such  vortex  induced  pressure  trough  has  been  discussed,  measured  and  computed 
previously  [2,10,13].  Figure  16  shows  the  comparison  of  the  measured  and  the  computed  static 
pressure  at  the  center  of  the  vortex  core.  The  computation  predicted  lower  values  than  were 
measured.  Figure  1 7  shows  the  trajectory  of  the  center  of  the  vortex  core  obtained  from  both 
computation  and  measurement.  The  computed  trajectory  agrees  with  the  measurement  very  well. 
Figure  18  shows  the  evolution  of  the  vortex  core  radius  obtained  from  both  computation  and 
measurement.  As  the  leakage  vortex  exits  the  cascade  passage  (X/ol.O),  the  shape  of  the  vortex 
core  becomes  asymmetric.  The  radius  of  the  core  is  an  averaged  value.  The  agreement  between 
the  computation  and  measurement  in  general  is  good.  The  measurement  data  shown  in  Fig.  16  to 
Fig.  1 8  are  extracted  from  reference  4. 

Figure  19  shows  the  comparison  of  the  computed  and  the  measured  static  pressure 
distribution  at  the  blade  surface  at  mid-span.  The  agreement  is  good  on  both  pressure  and  suction 
sides.  Figure  20  shows  the  static  pressure  distribution  near  the  blade  tip  ( 1 .5%  span  from  the  tip). 
On  the  pressure  side,  the  agreement  between  the  prediction  and  measurement  is  good,  however, 
there  is  a  pronounced  discrepancy  in  comparison  on  the  suction  side.  The  measurement  data 
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shown  in  Figures  19  and  20  are  extracted  from  reference  3.  Figure  21  shows  the  blade  to  blade 
profiles  of  axial  velocity  at  x/c=0. 1 1 .  The  agreement  between  the  prediction  and  the  measurement 
is  good  at  both  near  tip  and  mid-span  regions.  The  location  x/c=0. 1 1,  at  where  the  comparison  is 
made,  is  upstream  of  the  tip  leakage  vortex  roll-up  point.  Similar  velocity  profiles  are  shown  in 
Figs.22  and  23,  for  the  axial  locations  at  x/c=0.55  and  x/c=0.99,  respectively.  At  both  locations, 
the  calculations  over  predicted  the  magnitude  of  the  velocity  on  suction  side  near  the  tip,  where 
the  flow  is  dominated  by  the  passage  of  the  tip  leakage  vortex.  At  mid-span,  the  calculated  and 
the  measured  velocities  agree  well.  The  measurement  data  shown  in  Figures  21  to  23  are 
extracted  from  reference  10. 


DISCUSSION  AND  CONCLUSION 

The  tip  leakage  vortex  flow  is  a  challenging  subject  for  both  experiment  and 
computation.  The  experimental  data  from  Kang  and  Kirsch  [3,4,5]  provided  excellent  information 
for  a  numerical  simulation  of  a  vortex  flow.  Detail  flow  separations  around  the  tip  were  recorded 
and  illustrated.  The  measurements  were  conducted  with  five-hole  probe.  Although  the  diameter  of 
the  sensitive  part  of  the  probe  is  small  relative  to  the  blade  pitch  (1.4%);  however,  it  is  not 
insignificant  when  compared  with  the  diameter  of  the  vortex  core  (10.8%  at  x/c=0.8).  Limited 
information  regarding  to  the  experimental  uncertainty  can  be  found  in  reference  3. 

For  numerical  simulation,  three  different  grid  sizes  were  used  in  order  to  support  and 
resolve  the  large  gradients  of  the  static  pressure  and  the  velocity  components  within  the  small 
vortex  core.  The  grid  sizes  are  1.5  million,  2.5  million  and  3.0  million.  The  difference  between 
the  resulting  static  pressure  at  vortex  center  obtained  from  the  calculations  with  2.5  million  and 
3.0  million  grid  points  is  within  1%.  In  final  calculation,  grid  distribution  around  the  vortex  is  at 
least  or  better  than  the  20x20  (pitchwise  and  spanwise)  distribution.  The  numerical  differencing 
scheme  representing  the  RANS  equations  is  third-order  accurate  for  the  inviscid  flux  and  second- 
order  accurate  for  the  viscous  flux.  Some  common  and  more  noticeable  characteristics  of  the  tip 
leakage  vortex,  such  as  rolling-up  of  the  shear  layer  and  the  strengthening  and  weakening  of  the 
vorticity,  are  numerically  simulated.  The  predicted  trajectory  of  the  center  of  the  tip  leakage 
vortex  is  closer  to  the  blade  surface  than  it  was  measured.  The  calculated  radius  of  the  vortex  is 
somewhat  smaller  than  it  was  measured.  Comparing  with  the  measurement,  the  calculation  also 
predicted  lower  static  pressure  and  higher  velocity  in  the  vortex  core. 

Overall,  in  comparison  with  the  measurements,  the  numerical  simulation  over-predicted 
the  strength  of  the  tip  leakage  vortex.  Yet,  the  developing,  the  transporting  and  the  diffusing  of 
the  vortex  were  simulated  reasonably  well.  A  successful  computation  requires  a  high-order 
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accurate  differencing  scheme,  a  support  of  dense  grid  and  a  sophisticated  turbulence  model  to 
provide  correct  diffusion  at  the  vortex  core. 
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Figure  3:  Ink-trace  visualization  on  the  blade  tip  surface. 
(Kang  and  Kirsch  [4]) 


(Kang  and  Kirsch  [4]) 


Figure  5:  Computed  limited  skin  friction  lines  on  blade  tip  surface. 
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Figure  1:  Particle  tracing  of  vortices  around  tip  gap. 


Figure  2:  Schematic  of  the  vortex  structure  around  the  tip  gap. 
(Kang  and  Hirsch  [4]) 
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Figure  6:  Cross  flow  vectors  in  tip  gap  at  x/c  =  0.46. 
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Figure  10:  Cross  flow  vectors  an  static  pressure  contour  around  tip  leakage  vortex  at  x/c=1.125 
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Figure  1 1 :  Particle  tracing  of  the  vortices  around  the  tip  gap  with  vorticity  magnitude  contours  at 

various  traverse  planes. 
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Figure  14:  Surface  static  pressure  contours  at  blade  tip  and  suction  side. 


Figure  15:  Secondary  flow  structure. 
(Kang  and  Hirsch  [3]) 
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Figure  19:  Static  pressure  distribution  at  blade  surface  at  midspan. 
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